Toolbox for entanglement detection and fidelity estimation 
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The determination of the state fidelity and the detection of entanglement are fundamental prob- 
lems in quantum information experiments. We investigate how these goals can be achieved with a 
minimal effort. We show that the fidelity of GHZ and W states can be determined with an effort 
increasing only linearly with the number of qubits. We also present simple and robust methods for 
other states, such as cluster states and states in decoherence-free subspaces. 
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Due to recent advances in quantum control many ex- 
periments aim at the generation and manipulation of 
multipartite quantum states [H, 0j & 0] • I n such an exper- 
iment, one typically aims at the creation of some pure en- 
tangled multi-qubit state \ip), which can further be used 
for some quantum information processing task. Due to 
the unavoidable noise, however, the produced state will 
be some mixed state g cxp , which may significantly differ 
from the desired state \ip). 

First, one may be interested in the fidelity of the pro- 
duced state, that is, to what extent the desired state was 
prepared. This quantity is given by 



= (V'lPexpIV') 



(1) 



and should ideally equal to one. In practice, however, it 
is lower, nowadays experiments with four or more qubits 
achieve typical fidelities between 0.5 and 0.9. 

Second, one may ask whether the prepared state was 
indeed genuine multipartite entangled. This means that 
the iV-qubit state g cxp shows entanglement effects which 
cannot be produced by N — 1 qubits and all parties must 
have participated in the creation of the state g exp [f| . To 
verify this criterion of success, entanglement witnesses 
can be used [1, i, 0, i, ©, ES El El ■ These are observ- 
ables which have by construction a positive expectation 
value on all unentangled (biseparable) states, a measured 
negative expectation value signals the presence of entan- 
glement. 

A generic witness for a pure state \ip) is given by 



W = al- \ip)(ip\ 



(2) 



where a is the maximal fidelity of for bisepara- 
ble states, a quantity which can be directly computed 
[s|. By construction, this observable is positive on the 
biseparable states. Furthermore, from the expectation 
value Tr(g cxp W) we can determine the fidelity as — 
a—Tr(g cxp W), highlighting the close connection between 
fidelity determination and entanglement detection (l3j . 

For the implementation, two questions are relevant. 
First, since quantum state tomography requires an ex- 
ponentially increasing effort, the state g exp is often not 



completely known and one has to ask how many mea- 
surements are required for the evaluation of Tr(g cxp W). 
Typically, only local measurements are possible, and an 
observable like 



M = <7fc <g) a k ... <g) (Tfe 



(3) 



is called a local measurement setting [1, 0, Q • Note that 
measuring all the coincidence probabilities of the 2 N pos- 
sible outcomes of M. gives also information about observ- 
ables like 1 ® a k . . . <S> a k , etc. A second question concerns 
the robustness to noise of the witness W. As a simple 
model, one may consider the target state \ip) mixed with 
white noise, g(p) — p\ip)(ip\ + (1 — p)l/2 N and ask, how 
large p has to be, in order that g(p) is detected by W. 

It has been shown in Ref. @ that for important classes 
of states, namely GHZ states and cluster states, one can 
modify the witness in Eq. ^ such that it requires only 
two measurement settings and detects noisy states for 
p > 0.66 (GHZ states), and this independently of the 
number of qubits. This shows that entanglement detec- 
tion and fidelity estimation get not necessarily more dif- 
ficult when increasing the number of qubits. 

In this situation, three question are of interest. First, 
the question arises how a witness as in Eq. @ can be im- 
plemented by local measurements, and how many mea- 
surements are needed. Second, the exact payoff between 
the robustness to noise of a witness and the number of 
required measurement settings is of interest. Finally, the 
question arises whether there are constructions of wit- 
nesses beyond the projector-based witness of Eq. 

In this paper we address all three questions. We first 
show how the fidelity of GHZ and W states can be de- 
termined with N + 1 (resp. 2N — 1) local measurements. 
Then, extending the results of Ref. we derive wit- 
nesses for cluster states which are robust against noise 
and still require a small effort. Finally, we present a wit- 
ness construction for the four-qubit singlet state, which 
is simpler to measure than the construction in Eq. J2J), 
and, surprisingly, more robust against noise. 

Let us start by discussing the local decomposition of 
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the projector based witness for the case of GHZ states, 

\G N ) = -j=(\0)* N + \l) 9N ). (4) 

A witness for this state is given by W = 1/2 — \Gn)(Gn\ 
3, with the term \Gn)(Gn\ being decomposed as 
2|GAr>(Giv| = |0)(0|® JV + |1)(1|® JV + |0)(1|» JV + |1)(0|® 7V . 
The terms |0)(0| 8Ar and |1)(1|® JV can be directly mea- 
sured with the setting (cr z ) since they correspond to 
eigenvectors of it. Using the fact that |0)(1| = (a x + 
io~ y )/2 and |1)(0| = (a x — ia y )/2 we can write the re- 
maining term X = |0) <1 + \1){Q\® N as 

k N 

x = ^i E (-D fc/2 E(gK <S> *-> (b) 

k even tt i—1 i—k+1 

i.e. as an alternating sum of all products of the Pauli 
matrices o~ x and o~ y with an even number of o~ y . Here, 
denotes the sum over all permutations of the qubits, 
which yield different expressions. 

Now we choose the N measurement settings 

M k = [cos{^-)a x + sm(^-)a y ]® N , k = l,...,N; (6) 

which are just measurements in the x-y-plane of 
the Bloch sphere with different angles, and obtain 
yit^-i {—l) k M.k — NX as can be checked by direct cal- 
culation [14| . 

Observation 1. For the determination of the fidelity 
of the N-qubit GHZ state and evaluation on the witness 
W = -2/2 — \Gn)(Gn\ the N + 1 measurement settings 
(cz)®^ and Mk are sufficient. 

The required number of measurements can also be 
determined as follows: X in Eq. contains in each 
term only an even number of o~ y and X is symmetric 
under exchange of the qubits. The family of operators 
of this type is described by [N/2\ + 1 parameters [HI]. 
If we choose an arbitrary angle ctk € (0,tt/2) and define 
Mf = [cosK)a x ±sin(a fc )cTj / ]® Ar ,thenOT fe = Ml+M^ 
corresponds to two local measurements and also contains 
in each term only an even number of o~ y and is symmet- 
ric under exchange of the qubits. Therefore, by choosing 
L-/V/2J + 1 different atwe can express any observable of 
the former type as a linear combination of the VJlk ■ Not- 
ing that we may also choose a& = corresponding to 
(& X )® N (and, for N even, also = n/2 corresponding 
to (c y )® N ), we see that any observable which contains in 
each term (in the Pauli matrix representation) only an 
even number of a y (and a z elsewhere) and is symmetric 
under exchange of qubits can be measured by N local 
measurement settings. This argumentation also shows 
that the measurement directions are not unique. 

In order to see how this helps for other states besides 
GHZ states, let us consider W states [l6j], defined by 



where denotes again the symmetrized version. 

When decomposing |Wjv)(Wjv| into local measure- 
ments, the diagonal terms can again directly be measured 
with the (a z )® N setting. The off-diagonal terms, can, us- 
ing the fact that |0f)(10| + |10)(01| = (a x a x + cr y a y )/2 
and |0)(0| = (1 + <r z )/2, be written as 

2 1-JV 2 N 2 N 

y = ~w~ ( £ ® a * (gX 1+ ^)+E (g) °* (g)(l+^) 

7T i—1 i—3 7T i—1 % — 3 

Each of the two terms in this sum contains only an even 
number of o~ x (resp. o~ y ) and a different observable (l+cx z ) 
elsewhere. Hence, as discussed before, we can measure 
each term with ./V measurements settings of the type 
Aif = [(cos akO x j y ±sina/s(l + a z )] rg>N , which are effec- 
tively measurements in the s-z-plane (resp. y-z-plane) of 
the Bloch sphere. Using the fact that one of these set- 
tings can be chosen to be (1 + a z ) N , which is effectively 
a measurement of (a z ) N , we can summarize: 

Observation 2. The fidelity of an N-qubit W state 
can be determined by 2N — I local measurements. 

Consequently, also witnesses of the type W = al — 
|TT'v)OT"iv| (and also the more general witnesses in 
Ref. |2|) can be measured with 27V — I local measure- 
ments. Note that Observation 2 has general conse- 
quences in view of the results of Ref. [l2l ]. There, it was 
shown that any pure multipartite entangled state can be 
brought by local operations close to the W state. Con- 
sequently, a witness of the type W — al ~ \Wn)(Wn\, 
can, after reversing the appropriate local operations, de- 
tect any entangled pure state. Then, it is proved that 
the W state can be measured with N 2 — N + 1 settings, 
implying that any pure entangled state can be detected 
by N 2 — N + 1 measurements (although the robustness 
to noise may be small). Observation 2 shows that only 
2N — I measurements are already sufficient. 

The optimal decomposition for the three-qubit W state 
with five measurements was already given in Ref. [f| . For 
the four-qubit case, we obtain the decomposition 

\W 4 ){Wi\ = ^(-2^111cT z -4^11cT z OT z 

- 6]T ta z a z a z - %{a z ) m - 2(a x ) m - 2(a y ) m 

7T 

+ E [(i+^+^r 4 + (i+^-^H) (9) 

a= x ,y 

which requires the seven measurements of (a x )® A : 
(cr y ) m , (cr z ) m and (a z ± a a ) m with a = x,y. 

Let us discuss cases, where even with a simple de- 
composition the determination of the fidelity requires 
too much effort. This may happen in multi- photon 
experiments, where the measurement of a setting like 
(a z )® N may require a data collection time of several 
hours. In this case, instead of measuring the witness 
W = al — \ip) one can measure a witness like W' = 
al — + P where P is a positive operator. Since 
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P has only positive eigenvalues, (W) > implies that 
(W) > so W' is a valid witness, and allows to estimate 
the fidelity via > a — Tr(W g e x P )- The question is 
how to choose P in order require only few measurements, 
while still obtaining a good bound on the fidelity. 

Now we show how to derive a sequence of witnesses 
which are more and more robust to noise, while still be- 
ing simple to implement. The idea is to start from a 
known witness of the type W' = al — \ip) (ip\ + P and to 
systematically subtract terms from P. We demonstrate 
this for the six-qubit cluster state, the techniques can 
be straightforwardly adapted to other graph states. The 
six-qubit cluster state is given by 



\C 6 ) = i(|000000)- 



|111000> + |000111> — 1111111>) (10) 



This state is a graph state, corresponding after local rota- 
tions to the H-shaped graph [J] . A typical witness for this 
state would be W = 1/2 — |C 6 )(C 6 |, however, the meth- 
ods presented above yield a decomposition with 16 mea- 
surement settings. This witness detects a cluster state 
mixed with white noise for p > 31/63 w 0.492. 

In order to write down witnesses which require less 
effort, note that the cluster state can be described by its 
stabilizing operators [4j, |l7[ • These are 

g x = i<r z l(* s )® 3 ; g 2 = (^(l)® 4 ; g 3 = l(a^ 2 (l)® 3 ; 
54 = (a x )® 3 la z t; g 5 = (tf 3 (a z f 2 t; g 6 = (t) m (a z f 2 ; 

(11) 

and \Cq) is the unique state fulfilling gi\C§) = \C§). 

As the stabilizing operators describe the state, they 
can be used for the construction of witnesses. Indeed, as 
shown in Ref. [9| , a witness is given by 



w 1= 3i-2(n^p)-2(n 



(12) 



where the two groups of stabilizing operators correspond 
to the two-colourability of the H-shaped graph [10]. This 
witness requires only two measurement settings, namely 
(vz)® 3 (cr x ) m and (cr*)® 3 ^)® 3 . It detects states mixed 
with white noise for p > 0.7142. In order to show that W 
is a witness one can directly calculate that Wi — 2W > 0, 
i.e., Wi = 2W + P [|. 

To improve this witness, let us look at P = VVi — 2W > 
0. It is easy to see that P/2 is a projector onto a 49- 
dimensional subspace. The idea is to find a product basis 
such that many vectors of this basis lie in this subspace. 
Then, these vectors can be measured with one setting 
and they can be subtracted from the witness. In our 
case, the desired basis turns out to be the computational 
basis, i.e., the eigenvectors of the measurement (<j z )® 6 - 

Indeed, if we define for three qubits A = 1 — 
|000)(000| - |111)(111| the operator P - 2A ® A has no 
negative eigenvalues. Therefore, the observable 



is a valid witness, and requires three settings 
[(<J z )® 3 (a x )® 3 , {<J x ) m {a z ) m and (a z ) m ] for its mea- 
surement. It detects states mixed with white noise for 
p > 0.579, which is a reasonable improvement. 

As a sidestep, note the interesting fact that the wit- 
ness W2 may be improved by local filters. Such fil- 
ters are invertible operators Fi on each qubit, apply- 
ing then F = ® fe=1 -Ffe to a state, g 1— ► FgF\ keeps 
the entanglement (or separability) of the state. For the 
case of complete knowledge of g it was demonstrated in 
Ref. 0, that one can alternatively transform a witness 
asW^ F^WF, and optimize over F, in order to obtain 
a negative expectation value. Interestingly, if we restrict 
our attention to filters of the type Fi = ct;|0) (0| 
the witness i^W^-F 1 still needs the same three measure- 
ments as VV2, showing that filtering can also be useful for 
the case of incomplete information. 

In order to improve W2 further, we consider other 
eigenvectors of P, namely |(/) 1 / 2 ) = \vk) d> (| 001) ± 
|llO))/V2;|0 3 /4> = |« fc )®(|OlO)±|lOl))/V5and|0 B/6 ) = 
\v k ) ® (|100) ± |011))/V2, where \v k ) = |111) (\v k ) = 
1 000)) for k odd (even) and the upper signs hold for 

odd k. Then we can define V\ = J2i=i l<fo)(0i| = 
B ® C — |(1 — A) ® A, where A is the same as before, 
B = |111)(111| - |000)(000| and 
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2\/3 



{ 2-^J + 



So we arrive at the witness 



W 3 =W 2 - 2V X 



(14) 



(15) 



W 2 



Wi - 2A ® A. 



(13) 



which requires five measurement settings, namely the 
three settings for W2 and the two settings for measur- 
ing B ® C, and tolerates noise as long as p > 0.543. We 
can go further by considering six other vectors \<pk) for 
k = 7, 12, which arise form the 1^) by swapping the 
qubits 1,2,3 with the qubits 4,5,6. For them, one can 
define a V2 similar to V\, and one arrives at the witness: 

W 4 = W 2 - 2{V X + V-2) 

= VVi - 1®A-A® 1 - 2{B®C + C® B) (16) 

Note that for measurement of W4 the setting (<r z )® 6 is 
not required anymore, so in total only six measurement 
settings are needed. The witness detects states for p > 
0.5, so it is nearly as efficient as the original projector 
witness W = 1/2 — \Ce)(Ce\, however, VV4 requires a 
significantly smaller effort. It is this witness, which has 
also been used in a recent experiment 

There is another feature which makes the witness W4 
interesting from a general point of view. Namely, it can 
be written as 

m = ^-\Ce)(C 6 \ + \C 6 )(C 6 \ (17) 

where \C e ) = (-|000000) + |000111) + |111000) + 
1 111111) )/4 is an orthogonal cluster state in a different 
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basis. It is immediately clear that W 4 is a witness, which 
is of a similar strength as the witness W = 1/2 — | C 6 ) (Cg | 
since \Cq) and \C§) are orthogonal. Also, it becomes 
clear why the witness W 4 requires only a moderate ef- 
fort for its measurement: when comparing |C6)(Ce| and 
|Ce)(C76| + \Ce)(Ce\ in the standard basis, in the lat- 
ter some off-diagonal terms are canceled. As we know 
from the analysis of the GHZ and the W state, such off- 
diagonal terms are typically difficult to measure. 

Interestingly, other simple witnesses in the literature 
can also be viewed in this way. For instance, for the 
two-qubit singlet state |i/ ;± ) = ( 1 01} ± |10))/V2 a witness 
would be W = 1/2 — \if~)(ip~\ which can be simpli- 
fied to W = 1/2 - \i>-){ip-\ + IV'+XV'+I = (1 + ^® 

o~ x + a y ® cr t/)/2, which is known from Ref. [Io|. For 
the four-qubit cluster state \C 4 ) = (|0000) + |0011) + 
1 1100) — 1 1 1 1 1) ) /2 we can consider again the orthogonal 
state IC4) = (-|0000) + |0011) + |1100) + |llll))/2 and 
then the witness W = 1/2 - \C 4 )(C 4 \ + \C 4 )(C 4 \. This 
observable requires four measurements and has, in a dif- 
ferent form, already been derived in Ref. [Tll | . 

Observation 3. For cluster states, one can derive se- 
quences of witnesses, which require more and more mea- 
surements for an implementation, while becoming more 
and more robust to noise. Some of them can be improved 
by local filters, even if only few measurements have been 
done. Furthermore, a promising ansatz for a robust wit- 
ness for a state \tp) which requires only a moderate effort, 
is to consider W = al — \ip)(il)\ + \<j))(<f>\, where \<f>) is or- 
thogonal to 

Let us finally discuss an example how one can derive 
entanglement witnesses beyond the projector based wit- 
nesses. We consider states in decoherence free subspaces. 
These are pure states \ip), which are invariant under a si- 
multaneous unitary rotation on all qubits, i.e. we have 
(U)® N \ip) = e^\ip) (or, equivalently [|V>) (ip\, (U)® N ] = 
0). An example is the four-qubit singlet state [18| 

1*4) = -^(|0011) + |1100)-i(|01) + |10))®(|01) + |10))) 

(18) 

A natural witness for this state would be W = 3/41 — 
|\&4)(\&4|. It detects the state for p > 0.733 and requires 
the measurement of 15 local settings [3]. 



In order to write down a simpler witness, note that due 
to the unitary invariance of \^ 4 ) we have (o _ fc)® 4 |4 , 4) = 
1*4), for k — x,y, z. That is, observables of the type gk = 
(cfc)® 4 are stabilizing operators of l^t). Therefore, they 
may be a good tool to construct entanglement witnesses. 
It should be noted, however, that \^ 4 ) is not the only 
common eigenstate of the gk, the product of two singlet 
states, \<j>) = ® is a different one. 

In order to distinguish between the \^> 4 ) and the prod- 
uct of two-singlets, let us consider 

V(a k ) = (a k ) m + 3[K)® 2 (1)® 2 + (1)® 2 K)® 2 ] 

- I [l(o-fc) (82 l + lo- k lo- k + o-kto-ut + a k {tf 2 o- k ] 

(19) 

Here, the coefficients have been chosen in such a way, 
that for the four-qubit singlet state all terms have the 
expectation value one, e.g. 3{* 4 | [(cr fc ) <s>2 (l)® 2 ] |^ 4 ) = 1, 
etc. Then we consider Q — T>{a x ) + T>(a y ) + T>(a z ) 
For the singlet state we have ("I^l Ql^) = 21. For a 
biseparable state g^ s , direct numerical optimization gives 
Tr(g hs Q) < 14.35. Therefore, W = 14.35 ■ 1 - Q is 
a valid witness. This witness has interesting proper- 
ties: It requires only three measurement settings (namely 
(c/c)® 4 , k = x,y, z) and detects noisy states for p > 0.683. 
Therefore, despite of being simpler than the projector 
based witness, it is more robust to noise, and we have: 

Observation 4. For states in decoherence- free sub- 
spaces, entanglement witnesses can be constructed from 
their stabilizing operators g k = (erfc)®^. These witnesses 
are independent of the projector based witness. 

In conclusion, we have shown how to estimate the fi- 
delity and detect the entanglement for several important 
families of states. Our constructions demonstrate that 
one can typically solve these problems with a small ex- 
perimental effort. Therefore, our results will be useful for 
future experiments on multipartite entanglement. 
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